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The  observations  reported  in  this  paper  are  the  first  satellite- 
based  measurements  of  the  characteristic  Ka  emission  lines  of  nitrogen 
and  oxygen  from  the  earth's  sunlit  atmosphere.  In  addition  to  the 
atmospheric  X-ray  line  emissions,  lower  energy  X-ray  emission  was  also 
observed.  The  two  primary  mechanisms  for  the  production  of  X-rays 
from  the  sunlit  atmosphere  are  the  (coherent)  Thomson  scattering  of 
solar  X-rays  from  the  electrons  in  the  atomic  and  molecular  constitu¬ 
ents  of  the  atmosphere  and  the  absorption  of  incident  solar  X- 
radiation  followed  by  the  emission  of  characteristic  K  X-rays  of  ni¬ 
trogen  (0.392  keV),  oxygen  (0.525  keV),  and  argon  (2.94  keV).  Only  a 
few  observations  of  atmospheric  X-rays  produced  by  these  mechanisms 
have  been  reported  previously.  A  brief  summary  of  early  results  from 
the  experiment  herein  described  has  been  published  (RUGGE  et  al., 
1979). 

In  1967  HARRIES  and  FRANCEY  (1968)  observed  atmospheric  X-rays  in 
the  2-5  keV  region  during  two  rocket  flights.  The  X-rays  were  attri¬ 
buted  to  the  above  two  production  mechanisms,  primarily  on  the  basis 
of  observed  angular  distributions  (made  with  moderate  angular  resolu¬ 
tion:  10s  x  35°  FWHM)  which  matched  an  approximate  atmospheric  scat¬ 
tering  theory*  Since  these  rocket  experiments  had  effectively  no 
energy  resolution,  they  were  unable  to  determine  what,  if  any,  contri¬ 
bution  the  fluorescently  excited  Ar  X  line  made  to  the  flux. 


Also  in  1967,  GRADER  et  al.,  (1968)  observed  atmospheric  X-rays 
with  energies  above  0.55  keV  with  a  collimated  proportional  counter 
flown  on  a  rocket.  Collimation  was  *  30®  vertically  and  ±  5® 
horizontally.  The  authors  attributed  the  radiation  to  oxygen  and 
nitrogen  K  lines  fluorescently  excited  by  solar  X-rays,  based  on  the 
altitudes  of  observation  (peak  intensity  was  at  130  km)  and  the  pulse- 
height  distribution  of  the  measured  radiation.  Although  the  threshold 
of  the  pulse  height  analyzer  was  above  the  X  lines  of  0  and  N,  the 
counter  resolution  was  sufficiently  broad  that  the  pulse  height  dis¬ 
tribution  could  be  meaningfully  compared  to  that  which  would  be  ex¬ 
pected  from  the  line  radiation.  Moderate  agreement  was  found. 

Recently,  SEWARD  ^t_  al.  (1976)  have  observed  the  earth's  atmo¬ 
sphere  (in  the  1.5-9  keV  range)  during  solar  flares  using  a  propor¬ 
tional  counter  with  an  effective  field  of  view  of  5.2®  (FWHM)  on  the 
Ariel  V  satellite.  They  observed  the  Ar  K  fluorescence  line  as  well 
as  a  scattered  solar  component,  apparently  Including  solar  X-ray  line 
emission  near  6.5  keV.  Their  observed  Ar  K  X-ray  flux  was  a  factor  of 
about  100  below  that  predicted  by  using  an  approximate  calculation  by 
AIKIN  (1970). 

While  these  reports  have  clearly  established  the  existence  of 
atmospheric  X-rays  resulting  from  solar  X-rays  striking  the  earth's 
atmosphere,  they  do  not  provide  a  basis  for  a  complete  understanding 
of  the  phenomenon.  Indeed,  the  little  theoretical  work  that  has  been 
done  (AIKIN,  1970)  is  in  substantial  disagreement  with  the  one 
measurement  with  which  it  meaningfully  can  be  compared.  Our  intent  in 


this  paper  to  provide  sufficient  data  and  analysis  so  that  a 
comprehensive  understanding  of  these  "earth-albedo”  X-ray  phenomena 
can  be  obtained.  In  Section  II  the  exact  theory  is  developed. 
Section  III  describes  the  observations  made  and  the  instrumentation 
used.  Section  IV  presents  an  analysis  of  the  data,  and  Section  V 
summarizes  Che  major  points  of  the  paper. 

II.  THEORY  OF  EARTH-ALBEDO  X-RAYS 

HARRIES  and  FRANCEY  (1968)  presented  an  approximate  method  for 
calculating  the  X-ray  flux  from  a  sunlit,  planar,  single-constituent 
atmosphere.  Fluxes  for  both  scattered  and  Ar  K.  characteristic  radia¬ 
tion  were  computed.  When  the  approximation  is  extended  to  an  atmo¬ 
sphere  having  a  number  of  constituents,  such  as  the  earth's  atmo¬ 
sphere,  it  gives  satisfactory  results  in  many  cases.  We  will  use  the 
Harries  and  Francey  approximation  for  some  of  the  analysis  in  Section 
IV.  AIKIN  (1970)  derived  expected  atmospheric  X-ray  intensities  based 
on  a  presumably  more  complete  theory.  However,  there  was  substantial 
disagreement  between  the  predictions  and  the  measurement  of  SEWARD  et 
al.  (1976).  Unfortunately  AIKIN's  (1970)  paper  does  not  present 
enough  detail  to  allow  a  check  of  his  calculations  to  be  made. 

A  complete  and  accurate  theory,  assuming  fluorescence  and  co¬ 
herent  Thomson  scattering  of  the  incident  solar  X-rays ,  is  conceptual¬ 
ly  straightforward.  Because,  in  the  relevant  energy  range,  the  ab¬ 
sorption  cross-section  is  much  greater  than  the  scattering  cross- 
section,  only  single  scatterings  need  be  considered.  Consider  the 
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geometry  in  Figure  1.  The  vector  JL^  lies  in  the  yz  plane  and  extends 
from  the  scattering  point  P  to  a  detector  outside  the  atmosphere, 

making  an  angle  d  with  the  z  axis,  the  local  vertical.  Ji,  the  vector 
from  the  sun  to  P,  makes  an  angle  $  with  the  z  axis,  as  shown.  It 

does  not  necessarily  lie  in  the  yz  plane.  Let  the  detector  at 

distance  L  have  a  small  rectangular  field  of  view  of  a  radians 
parallel  to  the  projection  of  Jd  on  the  plane  normal  to  L^  (which  is 
also  the  plane  of  the  detector)  and  6  in  the  perpendicular 
direction.  The  box  at  P,  shown  in  the  expanded  view,  has  one  face 
normal  to  ^L  of  size  La  by  Lf3  (La  parallel  to  the  projection  of  H  on 
the  fac§  in  question).  The  third  dimension,  parallel  to  ^L,  is  dL 

(i.e.,  secQdz). 

A.  Scattered  X-Rays 

We  first  consider  the  scattered  flux.  If  the  flux  density  per 
unit  energy  of  solar  photons  incident  on  the  atmosphere  is  dNQ/dE 
(photon  -  cm-2  -s-1  -  keV“*),  then  the  flux  through  the  top  of  the  box 
(photon  -  s'1  -  keV_i)  is 

""df*^  -  d--~—  exp(-TI(z,E))(LS)(sec8dz)sinf,  (1) 

where  ¥  is  the  scattering  angle  (the  angle  between  Jd  and  L^,  secQdz  ■ 
dL,  and 


tI(z,E) 


(E)  / 


sec$ 'n^z'  )dz ' 


(2) 
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SURFACE  OF  EARTH 


Figure  1 


Geometry  for  a  satellite  observation  of  X-rays  from  the  sun-* 
lit  earth's  atmosphere*  The  interaction  region  is  expanded 
at  the  center.  is  the  vector  from  the  sun  to  the  interac¬ 
tion  point,  P,  which  is  at  an  altitude  z  above  the  earth's 
surface,  j*  is  the  vector  from  P  to  the  satellite,  and  the 
scattering  angle  ?  is  the  angle  between  Jg,  and  JL.  The  z  axis 
passes  through  P  and  the  center  of  the  earth.  The  other 
symbols  are  explained  in  the  text. 


The  angle  $'  is,  at  each  point  on  H,  analogous  to  at  P,  o.(E)  is  the 

**  i 

total  (absorption  plus  scattering)  X-ray  interaction  cross-section 
(HENKE  and  SCHATTENBERG,  1975)  for  atmospheric  constituent  i,  and  z  is 
the  altitude  at  which  the  interaction  occurs.  The  flux  density  at  the 
detector  (photon  -  cm  -s  -  keV  )  arising  from  scatterings  in  the 
box  is 


dNg(E) 

dF“ 


1  dF(z,E)V^/dq(E,f)s 
p-  dE  /  >  dfl  ;s,i 


-t0(z»E) 

ni(z)e  (Locscf), 


(3) 


J  **(  V  Of  A 

where  ( — — -)  is  the  differential  cross-section  for  coherent 

du  s ,  i 

scattering  (COMPTON  and  ALLISON,  1935),  and 

T0(z’E)  -L°i(E)  /  secfl’  ni(z')dz*  ,  (4) 

where  z^  is  the  altitude  of  the  detector.  In  equation  (3)  Lacsc’l'  is 
the  path  length  of  the  incident  photon  through  the  box.  The  total 
flux  density  is  obtained  by  integrating  equation  (3)  along  the  line  of 
sight: 


dNg(E)  dNQ(E) 


dE 


dE 


vec92I 


v  dfl  's  ,i 


(5) 


zr 

/  ni(z')  exp[-TI(z' ,E)  -  t0(z',E)J  dz’ 
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In  writing  equation  (5)  we  have  used  the  relationships  ft  ■  a&  and  dL  * 
secQdz.  Equation  (5)  holds  for  what  we  will  call  the  "penetrating" 
case,  the  case  wherein  a  backward  extension  of  L  strikes  the  earth. 
For  the  other,  the  "grazing"  case,  one  must  integrate  from  <*  to  zm  and 
then  from  zffl  to  z^,  where  zm  is  the  minimum  z  along  the  line  of  sight, 
j,.  The  optical  depths  in  equations  (2)  and  (4)  are  evaluated  by  using 
the  product  of  n^(z),  the  scale  height  of  constituent  i,  H^  and  the 
Chapman  integral  (e.g.,  SMITH  and  SMITH  1972). 

B.  Fluorescent  X-Rays 

Fluorescently  excited  X-rays  arise  following  absorption  of  inci¬ 
dent  X-rays  having  energy  above  E^  j,  the  K  absorption  edge  energy  of 
the  constituent,  j,  emitting  the  X-rays.  The  photon  flux  density, 
Nj.  j  (photon  -  cm-^  -  s-i)  is  derived  analogously  to  equation  (5): 


In  equation  (6)  <j  ( E )  is  the  cross-section  for  K  shell  photoelectric 

J 

absorption;  w  ,  the  X  fluorescence  yield  of  species  j,  is  the  frac- 
J 

tion  of  K  shell  absorptions  for  which  an  X-ray  rather  than  an  Auger 
electron  is  emitted  (TAWARA,  ^t  al. ,  1973);  and  Ep^j  Is  the  energy  of 
the  emitted  characteristic  X-ray. 
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Equations  (5)  and  (b)  are  not  amenable  to  analytic  solution.  In 
some  applications  to  be  discussed  below  the  expressions  are  too  un¬ 
wieldy  to  be  handled  by  machine  in  a  practical  manner.  Frequently  in 
the  grazing  incidence  case  the  approximation  ■  0  is  valid,  and 
the  equations  can  be  readily  calculated  by  computer.  For  the  pene¬ 
trating  case  we  have  extended  the  HARRIES  and  FRANCEY  (1968)  method 
for  use  with  a  multiple  constituent  atmosphere. 


III.  OBSERVATIONS 

The  observations  described  herein  were  made  with  the  two  low- 
energy  detectors  (LEU)  of  the  A-2  experiment  flown  on  the  HEAO-1 
spacecraft  launched  in  mid-August,  1977.  The  spacecraft  and 
instrumentation  wefe  designed  for  observations  of  cosmic  X-rays  and 
gamma  rays.  The  spacecraft  rotates  with  a  period  of  about  30  minutes 
about  the  sun-satellite  vector  H.  Since  L  is  normal  to  H  (i.e.  V  m 

A#  A» 

90°),  the  instruments  scan  the  earth  on  almost  every  rotation,  making 
possible  the  "accidental"  observations  described  herein.  The  LEDs  and 
the  HEAO  A-2  experiment  are  described  by  ROTHSCHILD  et  al.  (1979). 
Basically,  a  single  LED  is  a  multi-anode,  multi-layer,  thin  window, 
flow  proportional  counter  with  a  multigrid  collimator  immediately  in 
front  of  it.  The  collimator  is  designed  so  that,  depending  upon  the 
use  of  certain  sets  of  anodes,  the  detector  field  of  view  is  either 
2.95°  x  1.55°  or  2,55®  *  2.80®  (a  *  0).  In  the  smaller  f ield-of-view 
case  the  counter  area  is  194  cm^,  and  in  the  larger  FOV  case  it  is  212 
cm^.  The  energy  range  of  sensitivity  is  0.15  -  3.0  keV,  and  spectral 
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analysis  in  this  range  is  carried  out  by  using  a  32-channel  pulse- 
height  analyzer. 

The  32-channel  pulse-height  analyzer  used  with  the  LEDs  allows 
the  individual  channels  to  be  narrower  than  the  intrinsic  resolution 
of  the  detectors.  Hence  spectral  analysis  is  limited  only  by  the 
detector  resolution.  We  might  expect  to  derive  two  kinds  of  informa¬ 
tion  from  the  spectral  data.  First,  we  can  determine  the  strengths  of 
the  two  components  of  the  radiation  and  thus  the  "fluorescence-to- 
scattering  ratio".  (It  is  tempting  to  call  this  the  line-to-continuum 
ratio,  but  the  scattered  component  consists  largely  of  lines  emitted 
in  the  solar  corona.)  This  ratio,  as  a  function  of  geometry,  can  be 
used  as  a  check  of  the  theory  of  X-ray  emission  from  the  sunlit 
atmosphere.  Once  the  phenomena  are  understood,  similar  measurements 
can  serve  as  a  monitor  of  solar  activity,  in  particular  soft  X-ray 
emission,  for  satellites  (such  as  weather  satellites)  constrained  to 
view  the  earth.  Second,  the  ratio  of  N  to  0  characteristic  radiations 
might  be  used  as  a  monitor  of  atmospheric  composition.  As  will  be 
seen  below,  this  use  is  limited  because  the  measured  grazing  incidence 
spectra,  for  which  fluorescence  is  relatively  strongest  with  respect 
to  scattering,  are  heavily  dominated  by  0  emission,  since  atomic  oxy¬ 
gen  is  abundant  in  the  outer  atmosphere  and  the  detector  efficiency  is 
higher  at  0.525  keV  (0  Ko)  than  at  0.392  keV  (N  Ko).  Current  detec¬ 
tors  cannot  resolve  the  0  and  N  lines,  and  the  small  amount  of  N  emis¬ 
sion  present  has  little  effect  on  the  observed  spectrum. 
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The  presence  of  fluorescent  radiation  is  demonstrated  in  Figure 
2,  which  shows  the  spectrum  of  a  grazing  observation  (zm,  minimum  z 
along  L,  ■  279  km)  made  on  17  August  1977.  The  error  bars  show  the 
detector  counting  rates,  and  the  solid  curve  shows  the  best-fit  com¬ 
bined  fluorescence  and  scattering  spectrum.  The  statistic  for  the 
fit  was  7.97  for  9  degrees  of  freedom.  The  other  two  curves  show  the 
best  fit  fluorescent  radiation  spectrum  (N/0  count  ratio  was  a  fitting 
parameter)  and  scattering  spectrum  (solar  coronal  temperature  T  was  a 
parameter)  that  add  up  to  the  solid  line.  The  strong  peak  at  0.5  keV 
must  arise  from  fluorescence;  the  solar  spectrum  cannot  account  for 
it.  The  peak  in  the  scattered  spectrum  at  0.57  keV  arises  primarily 
from  strong  emission  by  the  0  VII  ion  in  the  corona,  with  a  contribu¬ 
tion  from  an  0  VIII  line  at  0.65  keV.  There  exists  some  uncertainty 
in  the  magnitude  of  this  emission  because  of  the  uncertainty  in  the 
solar  oxygen  abundance  and  the  high  temperature  sensitivity  of  the  0 
VII  emissivity  below  2  x  10^  K.  This  uncertainty  can  have  some  effect 
on  the  N/O  line  ratio  derived  from  a  spectrum,  but  the  solar  oxygen 
emission  fs  small  compared  to  the  atmospheric  fluorescence,  so  little 
error  in  the  fluorescence/scattering  ratio  is  introduced. 

The  measured  X-ray  intensity  as  a  function  of  the  sun-earth- 
satellite  geometry  gives  one  of  the  best  tests  of  the  albedo  X-ray 
theory  as  well  as  providing  a  potential  means  of  monitoring  solar 
activity.  The  accurate  HEAO-1  pointing  and  ephemeris  determinations 
allow  detailed  atmosphere  location  (latitude,  longitude,  altitude, 
6,  and  zm)  to  be  obtained  in  a  straightforward  manner. 
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Figure  2.  The  X-ray  spectrum  of  the  sunlit  earth  for  an  observation 
with  "grazing"  geometry.  The  best-fit  fluorescent  line  and 
scattering  components  are  shown  separately. 


In  Figure  3  the  (energy)  integrated  albedo  X-ray  intensity  (both 
scattered  and  fluorescent)  is  plotted  as  a  function  of  time  as  the 
satellite  scans  the  earth.  HEAO-1  was  at  local  noon  near  the  middle 
of  the  plotted  period  so  that  the  earth  was  viewed  almost  continuously 
during  the  30-minute  rotation  period.  The  flux  from  the  sunlit  earth 
was  about  two  orders  of  magnitude  above  the  diffuse  cosmic  X-ray  back¬ 
ground.  The  times  at  which  the  instrument  was  viewing  the  earth's 
atmosphere  rather  than  the  sky  are  indicated  in  the  figure  ("occulta- 
tion").  The  symmetrical  appearance  of  the  two  scans  results  from  the 
unusual  combination  of  the  phase  of  the  satellite  rotation,  the  ratio 
of  spin  period  and  orbital  period,  and  the  occurrence  of  local  noon  at 
T  •>  18  minutes. 


A  simplified  model  of  the  earth  albedo  X-ray  flux  is  helpful  in 
understanding  the  integrated  counting  rates  in  Figure  3.  We  have 
derived  equations  for  the  scattered  and  fluorescent  X-ray  flux  densi¬ 
ties  from  a  plane,  uniform,  optically  thick  atmosphere.  The  deriva¬ 
tions  are  similar  to  those  for  equations  (3)  and  (6)  and  will  not  be 
given  here.  In  many  cases  the  plain,  uniform,  optically  thick  atmo¬ 
sphere  is  a  useful  approximation  when  proper  account  is  taken  of  the 
composition  at  the  height  where  most  of  the  interactions  occur.  For 
scattering, 


dN  (E‘> 
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and  for  fluorescence, 
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Figure  3.  The  total  X-ray  counting  rate  plotted  as  a  function  of  time 
during  two  satellite  scans  across  the  sunlit  earth.  Local 
noon  is  at  T  ■  18  minutes  in  the  figure.  The  shaded  portions 
are  periods  during  which  spectra  were  accumulated. 
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The  behavior  of  n^  as  a 


where  Che  symbols  are  as  previously  defined 
function  of  z  can  be  approximated  by  the  form  nge-2^2^  where  H(z)  is 
the  (local)  scale  height.  For  this  reason,  at  least  in  the  penetrat¬ 
ing  line  of  sight  case,  the  bulk  of  the  X-rays  come  from  a  narrow 
height  range.  The  n^  in  equations  (7)  and  (8)  must  correspond  to  this 
interaction  height.  Note  that  the  fluxes  depend  upon  the  relative 
concentrations  of  the  atmospheric  constituents  rather  than  the  abso¬ 
lute  concentrations.  This  is  because  the  model  assumes  that  the  atmo¬ 
sphere  is  optically  thick.  HARRIES  and  FRANCEY  (1968)  presented  for¬ 
mulas  similar  to  equations  (7)  and  (8)  for  a  single  constituent  atmo¬ 
sphere. 


The  curve  in  Figure  3  can  be  understood  qualitatively  in  terms  of 
equations  (7)  and  (8).  During  the  rapid  rise  starting  at  I  *  2 
minutes,  9  ■  90°,  .but  the  optical  depth  of  the  atmosphere  along  the 
line  of  sight  is  not  large;  here  the  equations  are  invalid.  The  flux 
increase  is  due  to  the  increase  in  optical  depth.  The  peak  at  T  ■  4 
minutes  occurs  because  9  *  90°  and  the  atmosphere  is  optically  thick 
along  this  line  of  sight.  The  fluxes  then  decrease  as  the  earth  is 
scanned  and  9  decreases.  During  this  period  $  is  also  decreasing. 
Finally,  at  about  T  ■  13  minutes,  9  *  90°  and  the  atmosphere  is  opti¬ 
cally  thick  and  a  second  maximum  is  reached.  The  small  residual  above 
the  diffuse  background  near  T  ■  18  minutes  (local  noon)  arises  from 
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Che  atmosphere  at  and  above  che  satellite's  height  as  the  instrument 
views  in  a  direction  parallel  to  the  earth's  surface. 

Figure  4  presents  spectra  acquired  during  the  corresponding 
shaded  periods  in  Figure  3.  As  we  mentioned  above,  the  relative  im¬ 
portance  of  fluorescence  in  terms  of  counting  rates  is  higher  for  the 
grazing  observations  (9  *  90°).  (A  glance  at  Figure  2  will  show  that 
fluorescence  is  important  in  all  these  spectra.)  Equations  (7)  and 
(8)  are  useful  in  understanding  spectra  2-5.  The  solar  spectrum  is 
steep,  so  most  of  the  fluorescent  X-rays  arise  from  incident  photons 

having  energies  only  slightly  above  The  ratio  o.(EP  *)/o.(E) 

J  J  r » J  j 

is  therefore  small  and,  for  i  +  j,  o.(E  .  )/o  (E)  is  only  slightly 

l  r  ,  j  i 

larger  than  1.  This  means  that  the  fluorescent  component  is  more 
sensitive  to  cos0/cos<f>  than  is  the  scattered  component.  Therefore, 
as  9  approaches  90°,  che  fluorescent  radiation  increases  relative  to 
the  scattered  component.  This  is  apparent  in  spectrum  5.  In  the 
optically  thin  limit,  T ^  and  tq  approach  zero,  simplifying  equations 
(5)  and  (6)  somewhat.  The  relative  importance  of  the  two  processes 
depends  upon  the  cross-sections,  the  integrated  amount  of  material 
along  the  line  of  sight,  and  the  incident  spectrum.  Spectra  1  and  6 
clearly  show  that  N  and  0  fluorescence  in  the  sunlit  earth's  atmo¬ 
sphere  are  important  even  when  the  sun  is  not  very  active  (the  Ottawa 
2800-MHz  solar  flux  on  24  August  1977  was  only  84  solar  flux  units  — 
Solar  Geophysical  Data,  397,  part  I). 
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OBSERVED  SPECTRAL  DISTRIBUTION  (Counts) 


IV.  ANALYSIS 


We  have  carried  out  the  analysis  of  the  X-ray  data  as  follows.  A 
computer  program  was  used  to  obtain  fits  to  the  observed  spectra  with 
solar  coronal  temperature,  emission  measure  (see  below),  total  fluo¬ 
rescent  line  flux,  and  the  ratio  of  N  line  counts  to  0  line  counts 
(hereafter,  "N/0  ratio")  as  fitting  parameters.  The  X-ray  spectrum 
incident  upon  the  atmosphere  was  obtained  from  a  solar  X-ray  spectrum 
model  based  on  calculations  by  WALKER  (1973),  KATO  (1976),  and 
LOULERGUE  and  NUSSBAUMER  (1973).  The  solar  spectrum  was  then  used 
with  equations  (6)  and  (8)  and  the  CIRA  65  Atmospheric  Model  5  at 
12:00  hours  to  compute  the  fluorescence/scattering  and  N/0  ratios. 
The  computed  ratios  were  then  compared  to  the  results  from  the  best- 
fit  program.  A  further  check  was  obtained  by  comparing  the  parameters 
characterizing  the  solar  spectrum  derived  from  observations  closely 
spaced  in  time  but  having  very  different  geometries  and  hence  dif¬ 
ferent  observed  counting-rate  spectra. 

The  X-ray  spectrum  of  the  (optically  thin)  solar  corona  has  the 

form* 


$(E,T)  -  CN2dVf(E,T)  , 
e 


(9) 


*f  may  have  minor  Ng  dependence,  but  in  the  present  analysis  this  is 
unimportant. 
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where  Che  emission  arises  from  an  element  of  volume  dV  at  temperature 
T  and  having  electron  density  Ng.  In  general  the  corona  is  not  iso¬ 
thermal,  but  we  have  made  the  simplifying  assumption  that  it  is  in 
deriving  dNg(E)/dE: 


dNQ(E) 

dE 


[/  N*  dV]  f (E,T)  . 


(10) 


The  integral  in  brackets  is  the  emission  measure,  and  the  function 
f(E,T)  is  a  complex  combination  of  line  and  continuum  emissions.  In 
the  energy  range  below  a  few  hundred  eV,  where  the  scattering  compo¬ 
nent  is  dominant,  the  solar  spectrum  is  poorly  known.  Very  few  meas¬ 
urements  have  been  made,  and  the  uncertainties  in  the  line  emission 
are  particularly  large.  Our  computed  spectrum  in  this  region  is  based 
mostly  on  the  calculations  of  KATO  (1976). 

The  fitting  program  was  implemented  by  computing,  for  each  obser¬ 
vation,  the  "scattering  efficiency"  e  (see  equation  (5)): 

s 

es  *  sec9y^(dq-^ll)st  /  L  dz’ni(2>)  expt-x^z'  ,E)  -tq(z',E)].  (11) 
i  0 

(For  grazing  cases  the  Integration  limits  were  changed  in  accordance 

with  the  remark  following  equation  (5).)  The  scattering  part  of  the 

dNQ(E) 

counting  rate  spectrum  is  then  given  by  the  product  — — — , 

where  and  A  are  the  detector  efficiency  and  area  and  t  is  the  ob¬ 
serving  interval,  40.96  s.  Because  of  the  long  observing  interval  and 
the  satellite  spin  rate  of  12°/minute,  the  geometry  underwent  con- 
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siderable  change  while  each  spectrum  was  accumulated.  Therefore  it 

was  necessary  to  obtain  an  average  e  by  integrating  over  the  observa- 

s 

tion  period.  The  fitting  program  obtained  best- fit  emission  measures 

and  line  strengths  for  each  point  on  a  grid  of  temperatures  and  N/0 

2 

ratios.  From  these  the  best  fits  were  chosen  with  a  minimum  x  cri¬ 
terion. 

In  general  the  fits  obtained  were  satisfactory.  The  best  fit 

temperatures  were  almost  always  1-2  x  10^  K,  in  the  range  expected  for 

the  nonflaring  corona  (ACTON  et  al.,  1972;  MCKENZIE  et  al.,  1978). 

Emission  measures  were  5-12  x  10“^  cm”^  at  1  x  10^  K  and  2-9  x  10^ 

cm”^  at  1.5  x  10^  K.  While  somewhat  larger  than  the  1.5  x  10^  cm-^ 

obtained  by  ACTON  et  al.  (1972)  for  the  "general  corona,"  excluding 

active  regions  (at  1.45  x  10^  K),  these  values  are  not  unreasonable 

when  the  active  region  contributions  are  considered.  The 

fluorescence/scattering  and  N/0  ratios  will  be  discussed  below.  In 

2 

many  cases  acceptable  values  of  x  were  not  obtained.  In  general,  for 

the  grazing  observati'  ns  where  the  counting  rates  were  low  and  the 

2 

statistics  poor,  acceptable  values  of  x  were  obtained,  but  the  pene- 

2 

trating  cases,  with  good  statistics,  had  high  x  •  This  indicates  that 
systematic  errors  were  dominant.  Possible  sources  of  these  errors 
would  be  lack  of  knowledge  of  the  solar  spectrum,  failure  of  the  iso¬ 
thermal  solar  corona  assumption,  and  errors  in  the  detector  threshold, 
efficiency,  or  resolution  calibration. 

The  fluorescence  contribution  was  calculated  by  using  equation 
(6)  or  (8).  In  the  grazing  cases  equations  (5)  and  (6)  had  to  be  used 
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and,  because  the  effective  height  observed  varied  substantially  during 
an  observation,  the  equations  had  to  be  integrated  over  the  changing 
field  of  view.  The  approximation  ”  tq  “  0  was  ®ade  to  keep  the 
integration  from  being  too  unwieldy.  This  is  a  good  approximation; 

was  small  in  all  cases  and,  in  the  worst  case  treated,  tq  did  not 
exceed  0.32  at  any  energy. 

For  each  analyzed  spectrum  the  coronal  emission  measure  as  a 
function  of  coronal  temperature  was  derived  by  comparing  the  calcu¬ 
lated  fluorescent  flux  with  the  observed  fluorescence  as  determined  by 
the  best-fit  program.  The  emission  measure  as  a  function  of  tempera¬ 
ture  was  also  determined  by  using  the  best-fit  scattering  spectrum  at 
each  temperature.  These  two  curves  are  plotted  in  Figure  5  for  four 
grazing  observations  on  17  August  1977.  The  derived  emission  measures 
are  strongly  dependent  on  temperature  when  fluorescence  is  used  for 
the  determination  and  weakly  dependent  when  scattering  is  used.  This 
is  because  f(E,T),  the  solar  spectrum,  contains  a  factor  e-2^^  which 
has  a  strong  temperature  dependence  for  E  >  0.399  keV,  the  minimum 
energy  X-ray  that  can  fluoresce  nitrogen  (kT  ■  0.129  keV  for  T  ■  1.5  x 
10^  X).  In  contrast  one  can  see  from  Figure  2  that  scattered  photons 
of  energy  below  0.3  keV  have  the  greatest  impact  on  the  spectral 
fits.  The  curves  in  Figure  S  are  closely  spaced  except  for  the  earli¬ 
est  “scattering"  curve.  From  the  intersections  of  the  curves  we  de¬ 
rive  a  coronal  temperature  of  1.1  x  lO*5  K  and  emission  measure  of  7  x 

10^  cm""^.  This  temperature  is  below  the  best-fit  temperature  in  most 

2 

cases  (see  Figure  2),  and  x  would  lie  in  the  range  2-3  per  degree  of 
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CORONAL  EMISSION  MEASURE  Icm  Jl 


17  AUGUST  1977 

- 17:49:41  UT 

- 18:16:18  UT 

-  20:47:10  UT 

- 21:19:15  UT 

•  SCATTERING 

□  FLUORESCENCE 


SOLAR  CORONAL  TEMPERATURE  I10u  K) 

2 

Figure  5.  Derived  solar  coronal  emission  measure  (/ N  dV)  plotted  as  a 
function  of  coronal  temperature*  These  t\?o  parameters  were 
derived  from  fits  to  the  calculated  scattering  and  fluores¬ 
cent  contributions  to  the  X-ray  spectrum  of  the  sunlit 
earth's  atmosphere* 


freedom.  There  is  no  apparent  reason  for  the  one  outlying  curve.  All 
of  the  observations  involved  scans  over  a  wide  altitude  range,  typi¬ 
cally  200-425  km.  Interestingly,  a  solar  subflare  that  peaked  during 
the  second  observation  was  reported  by  a  single  observatory  (Solar 
Geophysical  Data,  397,  part  I).  The  flare  had  no  soft  X-ray  emission 
that  could  be  detected  by  observing  the  earth's  atmosphere. 

The  N/O  ratios  determined  from  the  best-fit  program  and  from  the 
integration  of  equation  (6)  are  shown  in  Table  1.  These  are  counting- 
rate  ratios;  the  detector  efficiency  is  taken  into  account.  A  range 
of  fitted  ratios  is  listed  based  on  a  952  confidence  interval  derived 
by  the  method  of  CASH  (1976).  The  calculated  ratios  are  generally 
considerably  smaller  than  the  best-fit  ratios.  We  have  previously 
discussed  possible  reasons  for  this.  In  summary,  for  the  grazing 
observations,  combination  of  the  fitting  program  with  calculated 
fluorescent  fluxes  yielded  a  coronal  temperature  and  emission  measure 
in  the  expected  range  and  fluorescent/scattering  ratios  consistent 
with  observation.  The  calculated  N/0  ratios  were  smaller  than  the 
observed  ones.  This  discrepancy  probably  arises  from  systematic 
errors  that  affect  the  fitting  program. 

For  the  optically  thick  observations  the  height  at  which  most  of 
the  observed  X-rays  originated  or  were  scattered  did  not  change  much 

during  an  observation.  Therefore,  the  uniform  plane  atmosphere  model 

» 

(equation  (8))  was  used  to  compute  the  fluorescence  and  cg  was  not 
integrated  over  the  changing  field  of  view.  The  calculations  were 
compared  with  the  best-fit  program  results,  and  some  observations  for 
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Table  1.  N/0  counting-rate  ratios 


UT  M/0  Ratio 

17  Aug.  1977  T  -  1.0  x  106  K  T  -  1.5  x  106  K 

Fit  Calc.  Fit  Calc. 


17:49:41 

0.50 

0.15 

0.25-0.50 

0.08 

18:16:18 

— 

0.33 

0.00—0. 50 

0.18 

20:47:10 

0.25-1.25 

0.15 

0.25-1.00 

0.08 

21:19:15 


0.75-1.00 


0.24 


0.50-1.50 


0.13 


which  integrated  counting  rates,  but  not  spectra,  were  available  were 
analyzed. 

Spectra  3,  4  and  5  of  Figure  3  were  analyzed  using  the  uniform 
plane  atmosphere  approximation.  For  each  spectrum  the  N/0  ratio,  the 
emission  measure,  and  the  fluorescence/scattering  ratio  were  deter* 
mined  from  the  fitting  program  and  by  calculation  for  T  *  1.0  and  1.5 
x  10^  K.  These  results  are  plotted  in  Figure  6.  The  intersections  of 
Che  graphs  can  be  used  to  specify  the  plotted  quantities.  The  coronal 
temperature  and  emission  measure  should  be  independent  of  geometry, 
and  in  fact  the  plots  show  little  variation  in  either  quantity.  We 
found  T  ■  0.9  -  1.4x  10^  K  and  emission  measure  »  7.5  x  10^9  cm~^. 
In  addition  an  analysis  similar  to  those  in  Figure  5  was  done  on  spec¬ 
trum  1  from  Figure  3.  We  found  T  *  1.1  x  10^  K  and  emission  meas¬ 
ure  •  8  x  10*9  cm“^.  These  results  agree  very  well  with  the  plane 
atmosphere  analysis.  We  point  out  that  the  fluorescence  calculations 
depend  upon  the  fluorescence  yields  uiK,  and  that  these  are  poorly 
known  for  0  and  N.  For  example,  we  use  values  of  TAWARA  et  al. 
(1973),  -  .00473  and  «  -  .00645,  while  FINK  et  al.  (1966) 

give  w.,  »  .0015  and  *  .0022.  Figures  5  and  6  indicate  that  a 

IvN  KU 

change  in  w  would  shift  the  "fluorescence”  curves,  resulting  in  a 

■V 

change  in  temperature  but  little  change  in  emission  measure  at  the 
Intersections. 

We  have  total  counting  rates,  but  not  spectra,  for  a  number  of 
observations  on  17  August  1977,  around  the  time  of  the  data  in  Figure 
5.  If  temperature  is  fixed,  emission  measure  can  be  derived  from  the 
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FLU0RESCENT/SCATTER1NG 


Figure  6.  Nitrogen/oxygen  line  counting  rate  ratio,  coronal  emission 
measure,  and  fluorescent/scattering  ratio  plotted  as  a  func¬ 
tion  of  coronal  X-ray  temperature.  The  parameters  were  de¬ 
rived  from  fits  of  the  observed  X-ray  spectrum  to  calculated 
model  spectra. 
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observed  counting  rates  and  the  plane  atmosphere  calculations*  The 
results  for  such  a  calculation  are  shown  in  Table  2.  For  1  x  10*3  K 
the  emission  measure  is  (6.4  ±  0.4)  x  10^  cm”^,  and  for  1.5  x  10^  K 
it  is  (3.1  ±  0.4)  x  10^  cm~^.  In  Figure  5  we  found  T  *  1.1  x  10**  K, 

A  Q  .O 

and  emission  measure  *  7.5  x  10  cm  •  It  is  clear  that  the  derived 
emission  measure  for  the  plane  atmosphere  calculations  is  somewhat 
smaller  than  that  for  the  grazing  observations.  Nevertheless,  at  1.0 
x  106  K,  near  the  temperature ’ found  for  the  grazing  observations,  the 
emission  measures  derived  under  the  plane  atmosphere  approximation  are 
tightly  grouped  for  a  wide  variety  of  geometries.  We  conclude  that 
the  uniform  plane  atmosphere  approximation  is  adequate  for  predicting 
the  X-ray  flux  from  the  sunlit  earth,  provided  that  the  region  in  the 
field  of  view  is  optically  chick. 


V.  SUMMARY 

We  have  analyzed  a  number  of  soft  X-ray  spectra  of  the  sunlit 
earth  obtained  by  the  HEAO-1  A-2  experiment.  The  X-rays  arise  from 
coherent  scattering  of  incident  solar  photons  and  fluorescence  of 
atmospheric  oxygen  and  nitrogen.  The  observations  were  compared  with 
theoretical  spectra  discussed  in  Sections  II,  III  and  IV.  Observa¬ 
tions  closely  spaced  in  time  but  having  a  wide  range  of  viewing  geome¬ 
tries  and  counting  rates  were  found  to  be  consistent  with  similar 
Incident  solar  X-ray  spectra.  The  derived  coronal  temperatures  and 
emission  measures,  parameters  characterizing  the  solar  spectrum,  were 
in  the  range  expected  for  the  nonflaring  sun  (ACTON  et  al.,  1972; 
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Table 

2.  Emission 

measures: 

plane  atmosphere  model 

UT 

6 

9 

Ealssloa 

Measure 

17  Aug.  1977 

(degrees ) 

(degrees) 

1.0  x  106  K 

1.5  x  106 

17:53:46 

84.6 

28.4 

5.4  (49) 

2.3  (49) 

17:54:27 

82.9 

27.3 

5.8  (49) 

2.5  (49) 

17:55:08 

80.8 

26.6 

6.1  (49) 

2.7  (49) 

17:55:49 

79.0 

26.4 

6.2  (49) 

2.6  (49) 

17:56:30 

77.6 

26.2 

6.1  (49) 

2.7  (49) 

17:57:11 

76.0 

26.2 

6.1  (49) 

2.7  (49) 

17:57:52 

74.0 

26.5 

6.1  (49) 

2.8  (49) 

18:09:53 

42.8 

54.3 

6.4  (49) 

3.2  (49) 

18:10:09 

42.9 

55.3 

6.4  (49) 

3.2  (49) 

18:10:50 

44.6 

57.8 

6.7  (49) 

3.3  (49) 

18:11:31 

47.1 

60.4 

6.6  (49) 

3.3  (49) 

18:12:12 

50.9 

62.8 

6.7  (49) 

3.3  (49) 

18:12:53 

55.4 

64.9 

6.7  (49) 

3.4  (49) 

18:13:34 

60.6 

67.2 

6.7  (49) 

3.3  (49) 

18:14:15 

66.2 

69.6 

7.1  (49) 

3.5  (49) 

18:14:56 

73.3 

71.8 

7.1  (49) 

3.4  (49) 

18:15:37 

83.9 

74.3 

6.2  (49) 

2.8  (49) 
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MCKENZIE  et  al.,  1978).  Therefore  It  appears  that  the  sunlit  earth's 
atmosphere  can  be  used  as  a  monitor  of  solar  activity  for  satellites 
unable  to  view  the  sun  directly. 

All  of  the  atmospheric  spectra  have  a  strong  contribution  from 
fluorescence.  Because  of  this  one  might  hope  to  monitor  atmospheric 
composition  by  detecting  characteristic  X-rays,  but  this  is  apparently 
not  practical.  A  major  obstacle  is  the  inability  of  proportional 
counters  (which  have  the  required  large  area  and  field  of  view)  to 
resolve  the  0  and  N  Ka  lines.  We  have  considered  the  use  of  gas -scin¬ 
tillation  proportional  counters,  which  have  resolution  approximately 
twice  as  good  as  that  of  a  conventional  proportional  counter  (e.g., 
ANDERSON  et_  al. ,  1978),  but  for  this  purpose  even  this  improved 
resolution  is  inadequate.  Additional  difficulties  arise  from  strong 
oxygen  line  emission  in  the  solar  corona  at  energies  near  that  of  0 
Ka  and  from  the  dominance  of  atomic  oxygen  high  in  the  earth's 
atmosphere  where  fluorescence  is  most  evident  (e.g..  Fig.  4,  Spectra  1 
and  6).  Remote  sensing  of  fluorescent  X-rays  can  be  a  useful  monitor 
of  atmospheric  constituents  only  Insofar  as  the  total  fluorescent  line 
counting  rate  (which  jLs_  predictable)  depends  upon  composition. 

X-rays  from  the  sunlit  earth  constitute  a  background  component 
that  might  interfere  with,  for  example,  satellite-based  auroral  X-ray 
observations.  Indeed  the  c}ata  reported  here  were  limited  by  the  deci¬ 
sion  to  turn  the  HEAO  A-2  experiment  off  when  viewing  the  earth,  in 
order  to  avoid  unnecessary  detector  fatigue.  In  a  separate  study,  we 
have  computed  the  X-ray  spectrum  from  the  sunlit  earth  during  a  period 


of  relatively  high  solar  activity  (but  in  the  absence  of  flares),  .  i 
compared  it  to  measured  auroral  spectra  (MIZERA  et  al.,  2978).  We 
found  that  even  weak  auroral  emission  should  exceed  the  solar  albedo 
flux  at  energies  above  about  2  keV,  even  with  unfavorable  geometry. 
At  energies  below  1  keV  the  solar  albedo  X-rays  may  cause  difficulty, 
particularly  for  observations  of  weak  aurorae.  Little  is  known  about 
the  auroral  X-ray  spectrum  at  such  low  energies. 


The  A-2  experiment  on  HEAO-1  is  a  collaboration  lead  by  E.  Boldt 
(Goddard  Space  Flight  Center)  and  G.  Garmire  (Cal  Tech)  with  collabo¬ 
rators  at  Goddard  Space  Flight  Center,  Cal  Tech,  Univ.  of  California, 
Berkeley  and  the  Jet  Propulsion  Laboratory.  This  work  was  supported 
at  The  Aerospace  Corporation  by  NASA  Contracts  NAS8-33235  and  NASW- 
3338,  by  the  Aerospace  Corporate  Research  Program  and  by  the  U.  S.  Air 
Force  Space  Division  Contract  F04701-79-C-0080  and  at  the  University 
of  California,  Berkeley,  by  NASA  Contracts  CIT  44-727600  and  NAS  5- 
23313.  We  thank  Mr.  Lee  A.  Christopher  for  developing  the  computer 
programs  used  in  the  analysis. 
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